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We study the performance of dynamical decoupling in suppressing decoherence caused by soft-cutoff Gaus- 
sian noise, using short-time expansion of the noise correlations and numerical optimization. For the noise 
correlations with odd power terms in the short-time expansion, there exists no dynamical decoupling scheme to 
eliminate the decoherence to arbitrary orders of the short time, regardless of the timing or pulse shaping of the 
control under the population conserving condition. We formulate the equations for optimizing pulse sequences 
for dynamical decoupling that minimizes decoherence up to the highest possible order of the short-time expan- 
sion. In particular, we show that the Carr-Purcell-Meiboom-Gill sequence is optimal in short-time limit for the 
noise correlations with a linear order term in the time expansion. 
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I. INTRODUCTION 

Quantum information processing [1] relies on the coher- 
ence of quantum systems. Unavoidable interactions between 
a quantum system and its environment (bath) introduce noise 
on the system and lead to error evolutions (decoherence) of 
the quantum system. Various methods have been proposed 
to combat the decoherence, including decoherence-free sub- 
spaces [2-4 1, error-correction codes [5 6], and dynamical de- 
coupling (DD) CHll . In particular, the DD scheme uses rapid 
unitary control pulses acting only on the systems to suppress 
the effects of the noise from the environments. DD has the ad- 
vantages of suppressing decoherence without measurement, 
feedback, or redundant encoding [8|. DD originated from the 
seminal spin echo experiment [10], in which the effect of a 
static random magnetic field (inhomogeneous broadening) is 
canceled. And more complex DD pulse sequences, such as the 
Carr-Purcell-Meiboom-Gill (CPMG) sequence Q31E], were 
designed to prolong the spin coherence time Ifl3l . 

The early DD schemes only eliminate low-order errors, 
i.e., the errors of quantum evolutions up to some low order 
in the Magnus expansion. By unitary symmetrization pro- 
cedure [7 8 1, DD cancels the first order (i.e., leading order) 
errors. To eliminate errors to the second order in short time, 
mirror-symmetric arrangement of two DD sequences can be 
used [8). The first explicit arbitrary Mth order DD scheme, 
which suppresses errors to 0(T M+l ) for short evolution time 
T, is the concatenated DD (CDD) JIIIIBI proposed by Khod- 
jasteh and Lidar. CDD sequences against pure dephasing were 
investigated for electron spin qubits in realistic solid-state sys- 
tems with nuclear spins as baths lfT6UT8l . Experiments |[T9T - 
23 1 have tested the performance of CDD. CDD works for 
generic quantum systems coupled to a finite bath Il24l l25l . 
However, since CDD uses recursively constructed pulse se- 
quences to suppress decoherence, the number of pulses in- 
creases exponentially with the decoupling order. As pulse 
errors are inevitably introduced in each control pulse in ex- 
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periments, finding efficient DD schemes with fewer control 
pulses is desirable. A remarkable advance is the Uhrig DD 
(UDD) Il26ti29l . UDD is optimal in the short-time limit in the 
sense that it suppresses the pure dephasing of a qubit coupled 
to a finite bath to the Mth order using only M qubit flips. The 
performance bounds for UDD against pure dephasing were 
established [ 30 1 . Shaped pulses OT1 [321 of finite amplitude 
can be incorporated into UDD Oil . Many recent experimen- 
tal studies [20, 23, 33-37 1 demonstrated the performance of 
UDD. 

It is important to find efficient schemes to suppress general 
decoherence (including pure dephasing and population relax- 
ation). Yang and Liu extended UDD to the suppression of 
population relaxation [ 29 1 . This inspired efficient ways to sup- 
press the general decoherence of single qubits, including con- 
catenation of UDD sequences (CUDD) [38 1 and a much more 
efficient one called quadratic DD (QDD) ||25l I39ti42l discov- 
ered by West et al 1391 . Based on the proof in (29), Mukhtar 
et al generalized UDD to protect arbitrary multilevel sys- 
tems with full prior knowledge of the initial states [43|. One 
can actually preserve the coherence of arbitrary multi-qubit 
systems by protecting a mutually orthogonal operation set 
(MOOS) l25l . By nesting UDD sequences for protecting the 
elements in the MOOS, the nested UDD (NUDD) Il25ll42ll43ll 
requires only a polynomially increasing number of pulses in 
the decoupling order. These universal DD schemes also work 
for analytically time-dependent baths l44l . 

The above-mentioned variations of UDD, however, rely on 
the finiteness of the baths, i.e., the existence of hard high- 
frequency cutoff in the noise spectra. Legitimate questions 
are: For quantum systems coupled to an infinite quantum bath 
or affected by soft-cutoff noise, can any DD be designed to 
eliminate the decoherence to arbitrary orders of precision in 
the short-time limit? And if yes, how can such DD be de- 
signed? Such questions have been previously addressed in 
some specific noise models. Comparing the efficiency of var- 
ious DD sequences in suppressing pure dephasing of a qubit 
due to classical noise, Cy winski et al observed that if the noise 
spectrum cutoff is not reached, CPMG sequences lfTTl[T2l ac- 
tually performs better than CDD and UDD sequences 11451 . 
Also, Pasini and Uhrig derived the equations for minimizing 
decoherence for power-law spectra, and found that the numer- 
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ically optimized sequences resemble CPMG [46]. Chen and 
Liu proved that for telegraphlike noise the CPMG sequences 
are the most efficient scheme in protecting the qubit coher- 
ence in the short-time limit and the decoherence can be sup- 
pressed at most to the third order of short evolution time by 
DD flTl . These results suggest that for noise with soft cutoff 
in the spectrum, there are certain constraints on the efficacy of 
DD. However, no conclusion has been drawn on the perfor- 
mance of DD with arbitrary timing and pulse shaping for the 
general cases of soft-cutoff noise. 

In this paper, we address the general question of the per- 
formance of DD against soft-cutoff noise, based on the short- 
time expansion of the noise correlations. The essence of time- 
expansion method is to avoid divergence in the integration of 
soft-cutoff noise in the frequency domain. We find that when 
the expansion of the noise correlation has odd-power terms 
in time, there exists no DD scheme to eliminate the decoher- 
ence to an arbitrary order of the short time, regardless of the 
timing and shaping of the control pulses under the popula- 
tion conserving condition. The odd-power terms in the time- 
expansion indicate that the noise has no hard high-frequency 
cutoff. Since for soft-cutoff noise the decoherence cannot be 
eliminated at a certain order of short time, we derive a set of 
equations to minimize the leading-order term in the short-time 
expansion while eliminating the lower orders. These equa- 
tions are numerically solved for optimal solutions. In particu- 
lar, for noise correlations with a linear order term in time, we 
prove that the CPMG sequences are optimal. For other noise 
correlations with odd-order terms, the minimum pulse interval 
of the optimized sequences is larger than in UDD sequences. 
This feature is important in realistic experiments when there 
is a minimum pulse switching time [48]. 

This paper is organized as follows: In Sec. [TTJ we analyze 
the performance of DD via the short-time expansion of corre- 
lation functions, and we give the condition under which deco- 
herence suppression to arbitrary order is impossible. The rela- 
tion between the high-frequency cutoff and the short-time ex- 



pansion of correlations is also discussed. In Sec. Ill we derive 
the equations for sequence optimization and obtain optimal 
DD for noise correlations with odd-power expansion terms. 



Finally, the conclusions are drawn in Sec. IV 



II. 



NO-GO THEOREM ON DYNAMICAL DECOUPLING 
AGAINST NOISE WITH SOFT CUTOFF 



We consider the pure-dephasing Hamiltonian for a single 
spin (qubit) 



(1) 



where cr z = |+)(+| -|-)(-| is the Pauli operator of the qubit, co a 
is the energy splitting of the qubit, and f3(t) describes random 
noise with average /?(?) = 0. Here the over bar denotes aver- 
aging over the noise realizations. We assume that the statistics 
of the noise fluctuations are Gaussian. 

After a duration of free evolution time T, the noise induces 
between the two states |±) a random phase shift J pit) that 



destroys the quantum coherence. We can suppress the deco- 
herence by DD control on the qubit. There is only one noise 
source /3(t) in the model Eq. ([T]), and to suppress the decoher- 
ence we need to protect a MOOS which consists of a Pauli 
operator cr x (more generally cr x cos (f> + cr y sin (f> with (f> being 
real) (25). We will prove later that DD can suppress the de- 
coherence (i.e., the protection of the MOOS \cr x }) only to a 
certain order of short evolution time for noise correlations that 
have odd-power expansion terms in time. We expect that the 
proof also applies to other quantum systems (e.g., multi-qubit 
systems) when the noise correlations have odd-power terms, 
since in those systems there are more noise sources and more 
system operators (e.g., a MOOS consisting of L > 1 Pauli 
operators) should be protected. 

When we apply a sequence of instantaneous unitary opera- 
tions cr x at the moments T\, Tj, ■ . ., T N , the controlled evolu- 
tion operator reads 



U(T) = {cT x ) N U{T N+u T N )(r x U(T N ,T N ^)--- 

where T$ = 0, T^+i = T, and the free evolution operator 

U(Tj +uTj ) = e-^£ +, ^ + ^ d '. 



(2) 



(3) 



Note that when N is odd, we may apply an additional cr x pulse 
at the end of the sequence for the identity evolution. Using 

r ,, T ^ -i%£ 0* 1 \ r a> a -m]dt ... 

o- x U(Tj +u Tj)o- x = e , (4) 

we write the evolution operator as 

U(T) = e 1 ^ £ ^FMdt e -if gp®FM* t (5) 

where we have defined the modulation function for instanta- 
neous 7r-pulse sequences 



FJt) 



(-iy 
o 



for? e(Tj,T j+l ] 
for t > T, or t < 



(6) 



Under DD control, the qubit is flipped at different moments, 
and the random field /3(t) is modulated by the modulation 
function F n (t), For multilevel systems, the modulation func- 
tions resulting from instantaneous 7r-pulse sequences may 
have values not restricted to {+1} for t e (0, T]. In Ref. [49], it 
is shown that for DD composed of specially engineered finite- 
duration pulses, the effective modulation functions can take 
values from [+1,-1,0) alternatively. We may also encounter 
effective modulation functions which are triangle wave func- 
tions during the time of system evolution ll50"l . For a more 
general analysis, we assume that the control conserves the 
populations and the phase modulation function F n {t) has a 
general form as 



F(t) = 



a bounded function for t e (0, T] 
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for t > T or t < 



(7) 
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At the moment T, the off-diagonal density matrix element 
of an ensemble is 



p u (T) = p n (0)e- ! '/o «.FW e -ii MFW t 



(8) 



where F(t) = F n {i) for ideal instantaneous <x v pulses on 
the qubit. The coherence is characterized by the ensemble- 
averaged phase factor 



W(T) = e-'h « f > F »<*. 



(9) 



For Gaussian noise, the ensemble-averaged phase factor W(T) 
is determined by the two-point correlation function P(t\)/3(t2) 
and W(T) becomes |45l l5fl |52) 



W(T) = e 



-X(T) 



where the phase correlation 



(10) 



(11) 



can be written as the overlap between the noise power spec- 
trum and a filter function determined by the Fourier transform 
of the modulation function P31 . The aim of DD design is to 
minimize x(T). 

We assume the noise is stationary, i.e., of time translation 
symmetry, f3(t\)f5(t2) - P{t\ - ?2)y6(0). The noise correlation 
function has the symmetry /3(t)fi(0) = /3(0)P(t), since /3(t) and 
/?(0) commute [For quantum noise, this may not be true. But in 
Eq. 
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, 1 1 and ?2 can be exchanged without changing the inte- 
gration. So the noise correlation can always be symmetrized]. 
These symmetries indicate that the noise correlation is an even 
function of time, i.e., 



p(t)J3(0) = C corr (f) = C corr (|f|). 



(12) 



The noise correlation can be transformed from the noise 
power spectrum S (w) as 



mm 



r 



da> . 
2^ 



S(a>)e 



(13) 



To analyze the decoherence suppression in short-time limit, 
we expand the noise correlation function in power series 



mm = YjC k 



(14) 



k=0 



The existence of odd-order terms means that the correla- 
tion function is non-analytical, which indicates that the noise 
source cannot be a finite quantum bath. The noise with non- 
analytical correlation functions must come from the fluctua- 
tions of an infinite bath, since otherwise the unitary evolu- 
tion of a finite quantum system will always lead to analyti- 
cal correlation functions. For example, the noise correlation 
/3{t)/3(0) = e~^' c has odd-order terms in the time expansion, 
and the noise has a Lorentz spectrum, which has a power- 
law decrease at high frequencies. We call this power-law de- 
crease a soft high-frequency cutoff in the noise spectrum 11461 . 



This kind of noise can be caused by Markovian (or instanta- 
neous) collisions in the bath ll53l . We can see from Eq. ( fT3| ) 
that when the noise spectrum has hard high-frequency cutoff, 
the noise correlation function /3(t)/3(0) is analytic as a func- 
tion of time and the expansion contains only even-order terms, 

AWcutoff = 2T=oC2*f 2 *. 

Expanding the noise correlation in time has the advantage 
of avoiding the divergence of integration over frequency for 
noise with a soft high-frequency cutoff in the power spectrum. 
Therefore we can derive certain general results on the efficacy 
of general dynamical decoupling control against soft-cutoff 
noise. Also, the performance of DD in the short-time limit is 
directly derived from the time-domain expansion. 

Using the expansion Eq. ( [14) , we write 



X(T) 



k=0 



(15) 



where the decoherence functions 



= [ dh f 

Jo Jo 



dt2{h-t 2 fF{t l )F{t 2 ). (16) 



Therefore if the modulation function F{t) is designed to make 
(0* = O}£ o = {0o = O,--- ,<f> M = 0}, (17) 

the decoherence in the short-time limit is eliminated to the 
Mth order, i.e., X (T) = 0(T M+3 ). Note that e~' £ = 1 

in Eq. (jSJ when 0o = 0. The even-order decoherence functions 
<f>2k describe the decoherence due to the low-frequency noise, 
while the odd-order functions 02/t+i arise from the soft high- 
frequency cutoff in the noise spectrum. 

A. Analysis of the decoherence functions 

Below we show that when the noise has soft high-frequency 
cutoff, it is not possible to suppress the decoherence to an ar- 
bitrary order, no matter how the DD modulation function F(t) 
is designed. 

The general filter function is defined as the Fourier trans- 
form of the general modulation function, 



F(t)e iwt dt = F*(-co). 

oo 

We expand it in power series 

(iaj) m - 1 



-A„ 



A m = m I F(t)t m - l dt. 



(18) 



(19a) 



(19b) 



Note that A m = 0(T m ). For a sequence of N instantaneous n 
pulses, F(t) is given by Eq. |6|l and A„, reads 

A£> = m J o F n (t)t m - X dt = £(-iy(7£ , - Tf ). (20) 
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As F (t) — for t ( (0, T], we extend the bounds of integra- 
tion for t to infinity and transform Eq. ( |T6] > to 



which is decomposed as (with the changes of summation order 
and indices) 



Pit = 



d(pi) 1 



xF(«!)F(w 2 )e 1 'e 22 <? 



(21) 



X°° /if,), f°° df,)n r°° T' 1 r*° J \2m+\2m-r+2 



where we set 77 — > after differentiation. Integrations over ?2, 
fi and g>i give 



0* = 



d(ir}} 



C doj F(to) - - 
* J_ M 27T ;(<w + if) 



Using 



A, + i, for r > 0, 



,=0 

and changing the summation index, we obtain 

,|2 k+1 



— 0. 



2^ 



(-ito) 



k+l 



-z 



(-z'w) 



i-r+2 „! 



(22) 



(23) 



(24) 



Note that the summation over r and the integration over fre- 
quency cannot be exchanged when the integration does not 



converge for each individual term. Using Eq. ( 19a 1, we have 



k\ 



f 



da> 
~2~n 



\F(u)\ 2 



k+\ k-r+2 

{-iio) k+l ~ 2j 2j 

r=l «=1 



(-l) k - KK 



(ico) 



k-r+3-n r \ n \ 



ifc+1 



z z 



{-\f- r KK 



k-r+3-n r \ n [ 



(25) 



r-l n=&-r+3 

We simplify the last line by using Eq. ( 19b| i and the equality 



n-r j 

t" = - — , for r > l,t > 0, (26) 



2;r ^ n\ 2 (r-l)! 

n-r ■ / 

which is proved in Appendix |A| We obtain 

\F(co)\ 2 ^^± 2 '-" k 



k\ 



f 



dot 
2^ 



(-ico) k+1 

k+l 



k+l k-r+2 . . . s _ r . . 

Zy (-1) A r A„ 



r=l 



r A r Ak-r+2 

r\ (jfc-r + 2)! 



(27) 



i. Even-order decoherence functions <p2 m 

We first consider the even-order decoherence functions 02m, 
which is caused by the noise with hard high-frequency cutoff. 
For even number 2m, , ^1; 
gral vanishes. Eq. (|27)i gives 



is an odd function and its inte- 



92;;/ 



£ Z Z 



2m+l 2m-r+2 



-(-iy A r A„ 

(/w) 2ffl - r+3 -" r! n! 



(2m) 



2m+\ 



- r A r A2 m - r +2 

r! (2w - r + 2)! ' 



(28) 



2/?7+l 2m-r+2 

z z 



r=l n=l 

-(-ir a,a„ 



-(-D/ A r A» 
2(/<y) 2m -'' +3 -" r! n! 



;-l n=l 
2m+l 



(2m)! 



z<-« 



2(;w) 2m -'' +3 -" r! n! 

2m-/- A2 m _ 



2m-r+2 



;-=l 



r! (2m -r + 2)! 



(29) 



The integrals of the terms with odd-power of co vanish. And 
for even functions of co, the sum n + r is an odd number, so 



(-l) r + (-l)" = 0. Thus the integral in Eq. (29 1 vanishes. With 



the change of summation index of the last line in Eq. (29 1, we 
write 



"2m 



= (2m)! 



(-D" 



A 



;;;+! 



(m+ 1)! 



7 -z<-» 



r A r A 



1m-r+2 



r! (2m -r + 2)! 

(30) 



From Eq. ( 30 1, we find that the following two sets of equations 
are equivalent 



[<hm = OCJo «* ( A m = 0}*!. 



(31) 



For instantaneous 7r-pulse sequences, the optimal solution 
of the equation set {A,„ = A„ = 0}^ =1 is 



Tsin 2 



2N + 2 



(j =1,2,..., JV), (32) 



which is the timing of UDD sequences B26I . Note that 
the decoherence functions 02m describe the noise with high- 
frequency cutoff in spectrum. Our result is consistent with 
Refs. 1 28 . 29]. The conditions Eq. ( flT) are more general than 
the one that leads to UDD, and Eq. \31\ may lead to more 
general designs of optimal DD. 



2. Leading odd-order decoherence function 1 



For the noise with soft high-frequency cutoff, the noise cor- 
relation function /3{t)/3(0) is not analytic. Its expansion con- 
tains odd-order terms C^+i^Zm+i + 0. Suppose the DD 
scheme has eliminated the leading order errors by making 

We want to study the leading 



! >2iu 



0} 



M-i 
;;;=() ' 



i.e., Eq. (31 



odd-order decoherence functions 4>2M-\ due to the soft high- 
frequency cutoff in the noise spectrum. 



The odd-order decoherence function <p2M-\ of Eq. (27i 
reads 



<p2M-\ 



(2M- 1)! 



2M 



z-< 



r A r A2Af+l-r 



(2M - 1) 

+ 2~K 

2M 2M-r+\ 

z z 



f 



da) 



r\ (2M+ 1 -r)\ 
\F{to)\ 2 



(-ia>) 2M 

i-iy a,.a„ 



{ico) 



2M-r+2-n r \ n \ 



(33) 



5 



The condition Eq. pi) gives 

4> 2M -i =(2M-1)! f 

J — i 

2M 2M+\-r 



2tt 



F(«)l 2 



(-/w) 2M 
(-l) r A r A„ 



(34) 



r=M+l n=M+\ 



(ioj) 



Notice in the summation 2M + 1 - r < M +1. We obtain 



(-1) M (2M- 



1) 



f 



da> \F(tJ) 
2n u> 2M 



(35) 



In Eq. (35 1, the integrand \F(ll>)\ 2 /a> 2M > and it cannot van- 
oo to oo. 



ish for all a> from 
theorem. 



Thus we have the following 



Theorem 1. For the noise correlation j3(t)/3(0) = zZT=o ^* l f l 
with \Cim + 0}^r ( J and C2M-I £ 0, there is no DD scheme that 
can induce a filter function F(t) to eliminate the errors up to 
0(T 2M+l ) for short evolution time T. That is, there is no DD 
scheme satisfying {4>2 m — 0125=0 and 4>im-i — 0. 

This theorem holds for arbitrary non-zero modulation func- 
tion F(t), and it shows that one cannot suppress the decoher- 
ence to arbitrary order when the noise correlation is not ana- 
lytic. For example, for the noise with the correlation function 
<?~ r ''', one cannot simultaneously eliminate the two leading de- 
coherence terms Q^o and C\(f)\, and the error induced by the 
noise is at least of the order 0{T y ). 

The result in Ref. B71 that DD can suppress decoherence at 
most to the third order of short evolution time for telegraphlike 
noise is general for noise with arbitrary statistics. In deriving 
Theorem[T] we have made the assumption that the statistics of 
the noise are Gaussian and the decoherence is determined by 
the two-point noise correlation. For non-Gaussian noise, there 
exists the decoherence effect from higher-order noise correla- 
tions. 



B. Noise correlation expansion and high-frequency cutoff 

Here we discuss further on the relation between the noise 
correlation expansion and high-frequency cutoff. Let us con- 
sider the correlation functions of the expansion form 



fc=0 



where the expansion coefficients 

1 d k C cm (0 



C k = 



k\ dt k 



are finite real numbers with 



C 2k -i = 0, for Jt < K, 
Ca-i * 0. 



(36) 



(37) 



(38a) 
(38b) 



The leading odd-order term in the short-time expansion of 
C CO rr(f) is C 2 K-i\t\ 2K ~ 1 ■ An example is C cou (t) = e" |f|3 with 



C\ — and C3 = -1. We assume that the noise correlation 
decreases smoothly at long correlation times, that is, 



lim — T C corr (f) = 0, for£ 
h» dt k 



0,1, 



and 



5 f 

Jo 



I L (co)= I e™— C COTr 



(t)dt, 



(39) 



(40) 



vanishes at large u> for L = 0, 1, . . . 11541 . The noise correla- 
tions that decay in the correlation time smoothly without fast 



oscillation satisfy Eq. ( 40 1. For example, the noise correlation 
e~ |f| has Il(d) = i(-l) L /(to + i) — * when a> — > 00. 

We consider the high frequency behavior of the Fourier 
transform of C corr (/), 



Xoo 
C cm ,(t)e^'dt. 

Integration by parts L > (2K +1) times gives 



S 2K (to) = 2% 



r=l v 



(r-l)l„ Ida)) 

Lj-i + 

(-icoY 



(-iu)) L 



(41) 



(42) 



where we have used Eqs. ( 37 1 and ( 39 1 



Using Eqs. (38 1 and (40 1, we obtain for large a>, 



Sik(<o) ~ 2 



(2K-iy. 

(iaj) 2K 



ClK-l + O 



1 



(43) 



which is a power-law decrease at high frequencies. 

When the noise correlation expansion contains only even- 
order terms, from Eq. (42 1 we have the noise spec- 
trum 5 even (w) = 2^R/l(<jj)/(-/w) l for an arbitrarily large 
L. From the assumption lim^co Il(oj) = 0, we have 
lim^oo 5 even (w)dLt L = for an arbitrarily large L and there- 
fore the noise spectrum has a hard high-frequency cutoff. One 
example is the correlation function e~'~ , which has the noise 
spectrum of exponential form ~ exp(-^-), and obviously the 
UDD sequence can suppress the noise effect order by order. 
The large a> behavior of other correlation functions of the form 
exp(- Yfj =y <X2jt 2 j) can be calculated by the saddle point inte- 
gration method, which gives a result of an exponential de- 
crease at high frequencies (i.e., hard cutoff). For example, 



when u> is very large, f°° e '' e lb "dt ^ \V> 

J & ' J-00 2 \ a{oj) 

g(co) = 3(f)5e- /27r/3 anda(w) = 12(f ) 2/3 e ! ' 2;r/3 . 
We now consider the noise spectrum 



Sp(o)) 



for to > Q, 



where 



(44) 



approximated by a power law (P > 0) decrease beyond the 
frequency Q. The noise correlation function is 



P(t)p{0) = C n ,cutoff(0 + C a ,p{t\ 
where the correlation due to high-frequency noise 

' da) 1 

Jn 



C n , P (t) = 2% f ——c 



(45) 



(46) 



6 



3' 


















o 





oo \ 



FIG. 1: The paths for the integral of l/z p 



and the noise correlation at low frequencies 
Cn,cutoff(0 = 251 



—S P (co)e-"' 

Jo 



2n 



(47) 



causes additional decoherence. The effect from Cn iCUt0 if(f) can 
be eliminated efficiently by DD sequences, such as UDD se- 
quences, as the noise spectrum of Cn iCUto ff(f) has a hard cutoff 
O. 

As Cn,p(f) is an even function, we just calculate the integral 
for the case of t > 0. We have (for t > 0) 



C Q At) 



7T 



1 



-,e m dz, 



where the paths cn, c,, and are shown in Fig. ([T]). 

Since the maximum of l/z p — > as \z\ — > °° in the upper 
half -plane, the contribution ^e m dz — 0. The contribution 



- f ±=j*dz=- r i l - p y- p e-y<dy 



(48) 



vanishes for even P. And 

r r 1 ,,. * 



^ L z pe w Jo nv™ 



e in,ei \i£l)e ie d6 



C£5 Cjj ' 



where 



/(1)= R y in^W/ rl ., | 
. r\ (r + 1 - p) ^ ' ' 

oo 



1! '^ki^Cr+l-p) 



1 / F f /> - 1 
Cn 2 (P - 1)! 



(49) 

(50a) 
(50b) 



r=0 



For even P, /£ n ' is an expansion of even powers of f, and 



,( 2 ) = 

e " 2(P-1)! 



for even P, 



(51) 



is the leading odd order expansion term of the noise correla- 
tion function. 

As an example, we consider the flowing noise spectrum, 



S 2 #( w ) = z ^ 

1 + L0 2K 



(52) 



for K e {1,2,...}, which can be transformed to a more gen- 
eral form of noise — -. % by using some scaling of the 
amplitude and a>. For example, the measured ambient noise 
for ions in a Penning trap has an approximate l/co 4 spec- 
trum at high frequencies and a flat dependence at low fre- 
quencies 11331 . which approximately corresponds to the noise 
spectrum Eq. ( |52| w ith K — 2. The corresponding correlation 
function of Eq. (|52|i is obtained by the inverse transform 



X°° da) e 
oo 2n~T+ 



Using the residue theorem, we have 



CO 



2K ' 



(53) 



mm 



2K 



i y exp[-/ e '^< 2 " +1 >|f|] 

exp[/^(2« + ' 



2K^exp[%(2n + l)(2K-l)] 
Expanding the right-hand side in powers of f, we get 

1 J^,{2n+\)(m-lK+\) 



2K 2K tt> - ! k 



(55) 



and summing the terms gives 



-i a {-i\t\r e ifc s^[(-ir + 1] 

m=0 



(56) 



In Eq. ( p6] l, the leading odd-order term of the time expansion 
is 2 (2g-i)! I ? l 2g-1 ' as predicted by Eq. ( jsij . By TheoremjlJ we 
have the following corollary. 

Corollary 1. Dynamical decoupling can only suppress the 
decoherence up to an error 0(T 2K+l ) in the time-expansion 
ofx(T)for the noise spectrum with high-frequency asymptote 
~ cj- 2k . 

For a more general noise spectrum with high-frequency 
asymptote to~ a , one can not suppress the error higher than 
the order of 0{T'" a+l ) with m a being the smallest even num- 
ber larger than a. This is because comparing with io~ m °, the 
asymptotic noise spectrum with a < m a decreases slower 
and the overlap between S(a>) and \F(a>)\ 2 is larger. 

The 1// noise is a special case (the case of P = 1) to 
be discussed. The noise correlation function Cn/> = i(f) = 
i cos z/zdz diverges at t — > 0; however, the decoherence 
function xp=i(T) [Eq. (JTTJ] induced by the high-frequency 
noise is finite. When there is no DD control, i.e., for the case 
of free induction decay, 



Xp=i(T) 



2nQ. 2 



[1 - cos(Qr) + Qr sin(Qr)] 



f 

Jm 



cosz 



dz, 



In Jn r z 
which has the short-time expansion 



(57) 



XP =l (T)=^(-2ln(Q.T) + 3-2y e )T 2 +0(Q. 2 T 4 ), (58) 
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with the Euler constant y e * 0.577. Notice the interesting 
result that^p = i(r) does not have the simple power law scaling 
at short time T. After eliminating the effect of low -frequency 
noise by high-order DD sequences which satisfy [A m = 0V™ =1 
(Eq. ( f3T| ) with 2M > P - 1, Xp(T) has a power law scaling 
^p(T) = 0(t p+1 ) EH, and we have Xp=i(T) = 0(t 2 ) when 
Ai = 0. 



III. SEQUENCE OPTIMIZATION IN SHORT-TIME LIMIT 

The aim of optimal DD is to minimize x(T). As indicated 
in Eq. ( 35 i, a smaller F(a>) at low frequencies will give a 
smaller 02M-1- Here we focus on DD with ideal instantaneous 
7i pulses. We use more pulses to construct a more efficient 
modulation function F(t) = F n {t) to minimize 02m-i, with 
the conditions {02,„ = OJ^Jo l^q. (31 1]. Using the method 
of Lagrange multipliers, we need to solve a set of nonlinear 
equations as 



bo- bo d cud abbot 1 aa 1 at! al □ d I b □ k> C 'ao d to hta b oi -and 



cd* □ a* no* cca aa ana 



oa □ OA 



c □ £> no jm io □ *o □ acq a>n 



AOa aq 



i ao at o □ ao c 



ODD 
UDD 
CPMG 



J I I I I 1_ 



J I 1_ 



_i I I I I L 



0.2 0.4 r £ l j r p 06 08 1 

FIG. 2: (Color online). Comparison of the ODD, UDD and CPMG 
sequences for different pulse number N. Squares (blue), triangles 
(green), and circles (red) correspond to UDD, CPMG, and ODD. The 
ODD sequences are optimized to minimize 3 under the constraints 
00 = 4>2 = 0. 



V{y,T}GM - 0, 
M 

G M = /^yj^f + 02M-1, 



/ 9 9 9 9 , 

VjyJ-j = ( , , , ). 

dTi dT N dyi dy M 



(59a) 
(59b) 

(59c) 



The introduced variables {yj} are the Lagrange multipliers. 
Note that the sequence optimization in Ref. [46 1 also used the 
constraints {A^' = 0), but the constraints were used there to 
guarantee the convergence of the calculation of x(T). Here, 
the constraints eliminate the lowest orders of errors ({02m = 
0}^Tq) [see Eq. ( [3l| ] in short-time limit. In particular, the 
decoherence from inhomogeneous broadening is eliminated 
when 0o = 0. 



We calculate Eq. ( 16 1 by separating the domain of integra- 
tion according to the value of F^t^F,,^)- For k > 0, we 
obtain 



-1 



(k + l)(k + 2) 
+ (T N+l 



N j-l 

4^^(r ; -r,)* +2 (-D 

>2 ,=1 



N 



N 

Then we have 



r,)* +2 (-i) 



N+l+i 



(60) 



dT k 



M 



2Yt(Tj-T k ) 2M (-iy-(T k -T ) 
j=k+\ 



2,1/ 



T k f M (-\) N+l 



k-i 

2YfTk-T]? M {-iy 



(61) 



For the special case of M = 1, G\ = y\ A ( ] ^+0 1 , we find that 
the CPMG sequences are solutions to Eq. ([59]). The timing of 



an -pulse CPMG sequence reads 
2j-l 



nCPMG 



2N 



T, for ; = l,...,N. 



(62) 



The CPMG sequences obviously satisfy the constraint A, 



[see Eq. (20 1], which is the so-called echo condition that 



eliminates the effect of static inhomogeneous broadening. 
Eqs. d6Tb and pOh give 



cm 
on 



(-D' 



k+l 



CPMG 



4N 2 



[l + (-!)% 



8T k 



2(-l) k+l yi . 



(63a) 
(63b) 



CPMG 



Thus the CPMG sequences also satisfy Eq. (59i with y\ = 

+ (~ l)j> so tne y are at ^ east tne l° ca lly optimal pulse 
sequences. It has been proved that CPMG sequences are the 
most efficient pulse sequences in protecting the qubit coher- 
ence against telegraph-like noise in the short-time limit l47l . 
With numerical evidence, we conjecture that it is also true that 
the CPMG sequences are optimal in the short-time limit when 
the time expansion of the noise correlation function has the 
two leading terms Co and Ci|?|. 

For other cases of minimizing 02m i with the condition 
{A^ = 0}^ =I , one can see that the short-time optimized DD 
(ODD) coincides with UDD for pulse number N < M. And 
as increases, the ODD sequences gradually approach the 
CPMG sequences. For example, ODD for the noise correla- 
tion 



mm = Co + c 2 t z + c 3 w + o(f\ 



(64) 



is shown in Fig.|2]in comparison with UDD and CPMG. The 
ODD sequences resemble the CPMG sequences when N is 
large. 

We show in Fig.[3ja) the performance of three DD schemes 
against the noise described by Eq. (64 1. A comparison is also 
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FIG. 3: (Color online). The decoherence function x(T) as a function 
of pulse number N (a) without hard cutoff and (b) with a hard high- 
frequency cutoff a) c = 40. Here the noise spectrum S(lj) = j^-t, 
with k = 10 5 and T = 0.5. Squares (blue), triangles (green), and 
circles (red) correspond to UDD, CPMG, and ODD. The ODD se- 
quences are the same as those shown in Fig. [2] 



shown in Fig.J^b) by considering a hard high-frequency cut- 
off u> c = 40. In Fig.[3ja), we can see that ODD sequences give 
better performance than UDD and CPMG sequences. These 
ODD sequences are optimal for a wide range of noise which 
has the noise correlation given by Eq. ( 64 1. When we intro- 



duce a high-frequency cutoff in the noise spectrum, as the case 
in Fig. [3jb), the ODD is the best initially when the number of 
pulses N < co c T/2 * 10, and the UDD sequences become 
better and suppress the decoherence order by order when N 
is large and the hard cutoff is reached. In Fig. [3jb), for large 
N UDD is better than ODD, since the ODD sequences are 
optimized for soft-cutoff noise rather than hard-cutoff noise. 
In Fig.|3ja) the decreasing of x(T) is a linear decrease in the 
double-logarithm plot, but in Fig. |3jb) it is much faster. This 
confirms that DD is not so efficient against soft-cutoff noise. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the dynamical decoupling control of deco- 
herence caused by Gaussian noise with soft cutoff. We have 
proved Theorem [T] which shows that, when the short-time ex- 
pansion of noise correlation has the (2K - l)th odd expansion 
term, DD can only suppress decoherence to 0{T 2K+l ). The 
existence of odd-order terms in the short-time expansion cor- 
responds to a soft high-frequency cutoff in the noise spectrum 
(i.e., a power-law asymptote at high frequencies). For these 
noise spectra, we have derived the equations for pulse se- 
quence optimization, which minimizes the leading odd-order 
decoherence function and eliminates even-order decoherence 
functions of lower orders. The ODD sequences obtained by 
this method coincide with the UDD sequences when the pulse 
number N is small, and they resemble CPMG sequences when 
N is large. For the special case that the short-time correlation 
function expansion has a linear term in time, the ODD se- 
quences are exactly the CPMG sequences. 

Although we derived Theorem [T] from a pure dephasing 



model, we expect that the result of the existence of the largest 
decoupling order in short-time limit can be extended to the 
general decoherence model (including both dephasing and re- 
laxation) of quantum systems. It is desirable to study the DD 
in suppressing the general decoherence of quantum systems 
in the soft-cutoff noise in the future. 
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Appendix A: Proof of Equation d26] 



To prove Eq. (26 1, we just need to prove 



lim f rfiV - — '—- — 

R ^J-R « ! 



(r-l)l 



forr>l, (Al) 



where the bounds in the integral guarantee that the modulation 
function F(t) is a real function. Using 



d y(±ix)_ = y R 

A n=r n=\ 

1 C R sinx ^ R" 1 

I ax — / 

— 1)!J_r x ^—^n\n(r — 



+ c.c), (A2a) 



0-1)! J_ R x "" ^ n\n(r-\)\ 

rR 

and litriR^oo J ^^dx = n, we just need to prove 
R" R" 



(i"- 1 + c.c), (A2b) 



lim y 



(n + r— l)ln n\n(r - 1)! 



(j"- 1 +c.c.) = 0. (A3) 



For r — 1, it obviously holds. For r > 2, we can show the 
difference 



A' (n 



R n (f-' + c.c.) 
(n + r — l)\n 



(r-l)\-Y](n + j) 



= <9 



a- 



(A4) 

so litriR^ooA = 0. By expanding the terms in the square brack- 



ets of Eq. ( A4 1, we have 



y *"(/"-' +C.C.) [ ^ ^ J 

Z-i (n + r-l)\ 

n=i v ' U=0 J 



a k n 



(A5) 



where is a number independent of n. We arrange the terms 
in the square brackets and get 



r-2 k 



k=\ 7=1 



(A6) 
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with bj independent of n. After some simplification it be- 
comes for r > 2 



Hence A = O(^), and Eq. (26 1 is proved 



(k+ifik Zj n \ 



/>=! 



+ c.c. = O 



(A7) 
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